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1. Introduction 



In hadron-hadron collisions, the most basic hard process is parton-parton scattering 
to form a large transverse momentum jet. The single jet inclusive transverse energy- 
distribution observed at the TEVATRON and CERN SppS shows good agreement 
with theoretical next-to-leading order O (a^) perturbative predictions over a wide 
range of jet transverse energies and tests the point-like nature of the partons down 
to distance scales of 10~ 17 m. However, data collected in Run I by the CDF collabora- 
tion at the TEVATRON indicated possible new physics at large transverse energy [1]. 
Data obtained by the DO collaboration [2] was more consistent with next-to-leading 
order expectations. However, because of both theoretical and experimental uncer- 
tainties no definite conclusion could be drawn. The experimental situation may be 
clarified in the forthcoming Run II starting in 2001 where increased statistics and 
improved detectors may lead to a reduction in both the statistical and systematic 
errors. 

The theoretical prediction may be improved by including the next-to-next-to- 
leading order perturbative predictions. This has the effect of (a) reducing the renor- 
malisation scale dependence and (b) improving the matching of the parton level 
theoretical jet algorithm with the hadron level experimental jet algorithm because 
the jet structure can be modeled by the presence of a third parton. Varying the 
renormalisation scale up and down by a factor of two about the jet transverse energy 
leads to a 20% (10%) renormalisation scale uncertainty at leading order (next-to- 
leading order) for jets with ~ 100 GeV. The improvement in accuracy expected 
at next-to-next-to-leading order can be estimated using the renormalisation group 
equations together with the known leading and next-to-leading order coefficients and 
is at the 1-2% level. 

The full next-to-next-to-leading order prediction requires a knowledge of the two- 
loop 2^2 matrix elements as well as the contributions from the one-loop 2-^3 
and tree-level 2 — > 4 processes. In the interesting large-transverse-energy region, 
E T ^> m quark , the quark masses may be safely neglected and we therefore focus on 
the scattering of massless partons. For processes involving up, down and strange 
quarks, which together with processes involving gluons form the bulk of the cross 
section, this is certainly a reliable approximation. The contribution involving charm 
and bottom quarks is only a small part of the total since the parton densities for 
finding charm and bottom quarks inside the proton are relatively suppressed. We 
note that the existing next-to-leading order programs [3, 4] used to compare directly 
with the experimental jet data [1, 2] are based on massless parton-parton scattering. 
Helicity amplitudes for the one-loop 2^3 parton sub-processes gg — > ggg, qq — > ggg, 
qq — > q'q'g, and processes related to these by crossing symmetry, have been computed 
in [5, 6, 7] respectively. The amplitudes for the six gluon gg — > gggg, four gluon-two 
quark qq — > gggg, two gluon-four quark qq — > q'q'gg and six quark qq — > q'q'q"q" 
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2^4 processes and the associated crossed processes computed at tree-level are also 
known and are available in [8, 9, 10, 11]. 

The calculation of the two-loop amplitudes for the massless 2-^2 scattering 
processes 

q + q^q' + q' (1.1) 

q + q->q + q, (1.2) 

q + q^g + g, (1.3) 

g + g^g + g, (1.4) 

has proved more intractable due mainly to the difficulty of evaluating the planar and 
non-planar double box graphs. Recently however, analytic expressions for these basic 
scalar integrals for massless particle scattering have been provided by Smirnov [12] 
and by Tausk [13] as series in e = (4 — D)/2. Associated tensor integrals have also 
been solved in [14] and [15] so that generic two-loop massless 2^2 processes can in 
principle be expressed in terms of a basis set of known two-loop integrals. With the 
notable exception of the maximal helicity violating two loop amplitude for gg — > gg 
which has recently been calculated by Bern, Dixon and Kosower [16] 1 , the two-loop 
matrix elements for the 2 — * 2 QCD parton scattering processes are not known. It 
is the purpose of this paper to provide dimensionally regularised and renormalised 
analytic expressions at the two-loop level for process (1.1) together with the time- 
reversed and crossed processes 

q + q' -»• q + q' , 
q + cf -> q + cf, 
q + q' -> q + cf. 

As is common in QCD calculations, we use the MS renormalisation scheme and 
conventional dimensional regularisation where all external particles are treated in D 
dimensions. We note that Bern, Dixon and Ghinculov [17] have recently completed 
the first full two-loop calculation of physical 2-^2 scattering amplitudes, the QED 
processes e + e~ — > fi + fi~ and e + e~ — > e~e + . There is an overlap between their QED 
calculation and the QCD results presented here and we expect that the analytic 
expressions presented here will therefore provide a useful check of some of their 
results. 

Our paper is structured as follows. In Section 2 we define our notation while a 
brief description of the methodology is given in Section 3. The results are collected in 
Section 4 where we provide analytic expressions for the interference of the two-loop 
and tree-level amplitudes as series expansions in e. Catani has described the pole 

1 This amplitude vanishes at tree level and does therefore not contribute to 2 — > 2 scattering at 
next-to-next-to-leading order O (af) . 
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structure of generic renormalised two-loop amplitudes [18] and we use his techniques 
to isolate the poles in the MS scheme. We find that the pole structure expected in 
the MS scheme on general grounds is indeed reproduced by direct evaluation of the 
Feynman diagrams. Ultimately these poles must be canceled by infrared singularities 
from tree level 2-^4 and one-loop 2^3 processes. The finite remainder of the 
two-loop graphs form the main results of our paper and are given in Section 4. Our 
findings are summarized in Section 5. 



2. Notation 



For calculational convenience, we treat all particles as incoming so that 

q(pi) + fa) + q'ips) + tip*) -> o (2.1) 

where the light-like momentum assignments are in parentheses and satisfy 

Pi + P2 + Ps + Pa = 0. 

As stated above, we work in conventional dimensional regularisation treating all 
external states in D dimensions. We renormalise in the MS scheme where the bare 
coupling ao is related to the running coupling a s = a s (fi 2 ) at renormalisation scale 
ji via 



ao S e = a 
In this expression 

S e = (47T) e e- e7 



2ttJ 



It 



OLs\ 

2n) 



+ O (afj 



7 = 0.5772 . . . = Euler constant 



(2.2) 



(2.3) 



is the typical phase-space volume factor in D = 4 — 2e dimensions, and (3 , Pi are the 
first two coefficients of the QCD beta function for N F (massless) quark flavours 

\1C\ - 10C a Tr N f - 6C F T R Np 



11^-47^^ 



6 6 
For an SU(N) gauge theory, where iV is the number of colours 

/iV 2 -l\ 1 
, C A = N, T R = -. 

2JV ' ' 2 



The renormalised four point amplitude in the MS scheme is thus 



\M) = 47ra s 



(2.4) 



(2.5) 



(2.6) 



where the \M^) represents a colour space vector describing the i-loop amplitude. 
The dependence on both renormalisation scale \i and renormalisation scheme is im- 
plicit. 
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We denote the squared amplitude summed over spins and colours by 

(M\M) = A(s,t,u), (2.7) 

where the Mandelstam variables are given by 

s = (Pi+P2)\ t = (p 2 +p 3 ) 2 , u =( Pl+ p 3 ) 2 . (2.8) 

For the physical processes, the spin and colour averaged amplitudes are related to A 
by 

! A(s,t,u) (2.9) 





>q> + q>)\ 2 


Y,\M{q + q> - 


^q + q')\ 2 


J2\M(q + q>- 


^q' + q)\ 2 


^2\M(q + q'~ 


^q + q')\ 2 



AN 2 

^A(u,t,s) (2.10) 

^A(t, 8 ,u) (2.11) 

^A(u,t,s). (2.12) 
The summed and squared amplitude has the perturbative expansion 

A\s, t, u) + (f) A %s, t, u) + (y-) 2 A\s, t,u) + (afj 



A(s,t,u) = 16n 2 a 2 
In terms of the amplitudes 



(2.13) 



A%s,t,u) = (M®\M<®) = 2(N 2 - 1) ( - e] , (2.14) 



A\s,t,u) = ((M {0) \M {1) ) + (M W \M (0) )) , (2.15) 
A 8 (s,t,u) = ((M {1) \M {1) ) + (M (0) \M {2) ) + (M {2) \M {0) )) . (2.16) 

Expressions for A 6 are given in Ref. [19] using dimensional regularisation to isolate 
the infrared and ultraviolet singularities. 

Here we concentrate on the next-to-next-to-leading order contribution A 8 and 
in particular the interference of the two-loop and tree graphs. 

3. Method 

Massless two-loop integrals for 2 — > 2 scattering can be described in terms of a basis 
set of scalar master integrals. The simpler massless master integrals comprise the 
trivial topologies of single scale integrals which can be written as products of Gamma 
functions: 



Sunset (s) = — ^ — (s) 
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Glass( S )= _OQ-( s ) 



Tri( S )= _(Q 



the less trivial non-planar triangle graph [20], 

Xtri(s) = 




and two scale integrals that are related to the one-loop box graphs [21, 22], 

TT(m) 



Abox(s,t) = 



Cbox(s,t) = /J (s,t). 
The planar double box [12] and non-planar double box [13] 



Pboxx(s, t) = 
Xboxx(s, t) = 



(s,t) 



involve multiple Mellin-Barnes integrals and are much more complicated to evaluate 
as series expansions in e. Expressions for these integrals valid through to O (e°) are 
given in [12] and [13] respectively. 

It turns out that for the two latter topologies, integrals involving loop momenta in 
the numerator cannot be entirely reduced in terms of the simpler integrals mentioned 
above and an additional master integral is required in each case. Reference [14] 
describes the procedure for reducing the tensor integrals down to a basis involving 
the planar box integral 



Pbox 2 (s, t) 



(s,t), 



where the blob on the middle propagator represents an additional power of that 
propagator, and provides a series expansion for Pbox 2 to O(e ). However, as was 
pointed out in [23], knowledge of Pbox! and Pbox 2 to O (e°) is not sufficient to 
determine all tensor loop integrals to the same order. A better basis involves the 
tensor integral, 

Pbox 3 (s,t) = [©1 ( (s,t), 



where © represents the planar box integral with one irreducible numerator associ- 
ated with the left loop. Symmetry of the integral ensures that, 



© 



[s,t) = 



© (s,t). 
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Series expansions for Pbox 3 are relatively compact and straightforward to obtain and 
are detailed in [24, 25]. Pbox2 can therefore be eliminated in favor of Pbox3. We 
note that this choice is not unique. Bern et al. [17] choose to use the Pboxi and 
Pbox 2 basis, but with the integrals evaluated in D = 6 — 2e dimensions where they 
are both infrared and ultraviolet finite. 

Similarly, the tensor reduction of the non-planar box integrals [15] also requires 
a second master integral, 



where the blob again denotes an additional power of the propagator. For the non- 
planar graphs there are no complications as in the planar case and all tensors to O (e°) 
may be described in terms of the series expansions of Xboxi and Xbox 2 through to 



In general tensor integrals are associated with scalar integrals in higher dimension 
and with higher powers of propagators. This connection can straightforwardly be 
achieved using the Schwinger parameter form of the integral and is detailed in [22] 
where explicit expressions for generic two-loop integrals with up to four powers of 
loop momenta in the numerator are given 2 . Systematic application of the integration- 
by-parts (IBP) identities [27] and Lorentz invariance (LI) identities [28] is sufficient 
to reduce these higher-dimension, higher-power integrals to master integrals in D = 
4 — 2e. Some topologies that occur in Feynman diagrams such as the pentabox [22] 
are immediately simplified using the IBP identities and collapse to combinations of 
master integrals. However, the tensor integrals directly associated with the master 
integrals usually require more care. Explicit identities relevant for the tensor integrals 
of the Abox and Cbox topologies are given in [22], for Pboxi and Pbox2 integrals 
in [14] while those for the Xtri, Xboxi and Xbox2 integrals are detailed in [15]. Using 
these identities, we have constructed MAPLE and FORM programs to rewrite two- 
loop tensor integrals for massless 2-^2 scattering directly in terms of the basis set 
of master integrals. 

The one-loop integrals are much better known. There are only two master inte- 
grals, the scalar bubble graph, 



2 A method to reduce tensor integrals constructing differential operators that change the powers 
of the propagators as well as the dimension of the integral was presented in Ref. [26] . 




O(e ) [13, 15]. 



Bub( S )= -Q-(s) 



and the one-loop scalar box graph, 



Box(s, t) 



(s,t). 
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We treat the tensor integrals in the same way as the two-loop integrals: shifting both 
dimension and powers of propagators and then using IBP to rewrite the integrals as 
combinations of Bub and Box. We note that this is not a unique choice for the 
master integrals. The one-loop bubble graph is proportional to the one-loop triangle 
graph with one off-shell leg. Another common choice is to replace the one-loop box 
in D = 4 — 2e by the finite one-loop box in D = 6 — 2e, Box 6 . 

The general procedure for computing the amplitudes is therefore as follows. First 
the two-loop Feynman diagrams are generated using QGRAF [29]. We then project by 
tree level, perform the summation over colours and spins and trace over the Dirac 
matrices in D dimensions using conventional dimensional regularisation. It is then 
straightforward to identify the scalar and tensor integrals present and replace them 
with combinations of master integrals using the tensor reduction of two-loop integrals 
described in [14, 15, 22] based on integration-by-parts [27] and Lorentz invariance [28] 
identities. The final result is a combination of master integrals in D = 4 — 2e which 
can be substituted for the expansions in e given in [12, 13, 14, 15, 21, 22, 24, 25]. 

4. Results 

In this section, we give explicit formulae for the e-expansion of the two-loop contribu- 
tion to the next-to-next-to-leading order term A 8 (s, t, u). To distinguish between the 
genuine two-loop contribution (M^\M^) + {M^\M^) and the squared one-loop 
part (M^\M^), we decompose A 8 as 



straightforward to obtain. For the remainder of this paper we concentrate on the 
technically more complicated two-loop contribution A s ( 2x0 ). 

We divide the two-loop contributions into two classes: those that multiply poles 
in the dimensional regularisation parameter e and those that are finite as e — > 



Voles contains both infrared singularities and ultraviolet divergences. The latter 
are removed by renormalisation, while the former must be analytically canceled by 
the infrared singularities occurring in radiative processes of the same order. The 
structure of these infrared divergences has been widely studied and, as has been 
demonstrated by Catani [18] , can be largely predicted. 

4.1 Infrared pole structure 

In the notation of Section 2, the universal infrared divergences present in a one-loop 
amplitude are given by the factorization formulae 




A 8 {2x0) (s, t, u) = Voles + Finite. 



(4.2) 



M {1) ) = I (1) (e) \M {0) ) + \M {1)&n ), 



(4.3) 



7 



where \A4^ &n ) is finite as e — > and the singular dependence is determined by the 
colour-charge operator /^(e) that acts on the tree- level colour vector For 
the n parton process we sum over all possible colour antennae with colour operators 
Tj • Tj acting on the state \M^) to obtain 



( ) 2r(l-e)^^ Tj ' Ve 2 + J V 2p,p, J (4 ' 4) 

and where A^ = — 1 if « and j are both incoming or outgoing partons and A^ = 
otherwise and the constants 7$ are given by 

7<? = 7? =2' 7 S = ^r- (4.5) 
Similarly at the two-loop level there is a factorisation of the infrared singularities 
\M {2) ) = / (1) (e) \M W ) + J (2) (e) |M(°)) + |7W^ fin > (4.6) 

where now 

I^(e) = - V(e) (lU(e) + + + k) J«(2 C ) + H«( C ) 

(4.7) 

where the constant X is 

'-(S-t)°-t^ <«» 

The function contains only single poles and is process dependent. For the case 
of the quark form factor (in the MS scheme) it is given by 

1 p e7 / 2 p -iAi 2 7r\ 2e 
JEf( 2 )( e ) = 7-—^ r ( ^ (4.9) 

1 ; 4er(l-e) I 2 Pl .p 2 i 1 ; 



with 

28 „ „ /16 



#(2) 



-7(i) + 3C F X + -C 2 A,C F - - (3 C F - (- - 7C 3 ) CjrC, 



(4.10) 



where Cn is the Riemann Zeta function, £2 = 7r 2 /6, £3 = 1.202056 . . . and 

7(i) = (-3 + 24C 2 - 48Cs) Cl + (-y - y( 2 + 24C 3 ) C F C A + (J + y( 2 ) C F T R N F . 

(4.11) 

We expect that in the four-quark two loop amplitude, we might obtain contributions 
from for each of the six colour antennae. 
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Applying the formalism to the case at hand, we find that the pole structure of 
the two-loop amplitude interfered with tree level has the following structure 



Voles = 2 Re 



l(M^\I W (e)I^(e)\M^} - ^ (M^\I W (e)\M^) 



+ (M {0) \I {1) (e)\M {1)fin ) 

- ^(1-26) 



+e 



+ (^(°)lif( 2 )( e )l^(°)) 

The colour algebra is straightforward and we find 
(M^\I^(e)\M w ) = (M^\M W ) 



^j + Kj (M^\I {1 \2e)\M^) 



on 



x 



1 3 
r(l-e) Ke 2 + Ye 



1 

A 



.EL 

s 



2 
A 



.EL 

u 



A 2 



A 



.EL 
t 



(4.12) 



(4.13) 



(M^\I W (e)I^(e)\M^} = (.M (0) |.M (0) > 

A 4 - 3 A 2 + 3 / 



x 



e 2e7 (\ 3\ 2 
r(l-e) 2 [e 2 + Ye) 
N 2 — 2 ' - 2 

~~ 2 A 2 \ .s 



AT2 

,2\ ' 



^\ 2e jy2 + 3 / ^2 
' t J + iV 2 I u 



2c 



4 

A 2 



s 



u 



+ 



N 2 



.EL 

s 



2c 



(4.14) 



and 



(M^\I W (e)\M W 



fin\ 



3 £7 



X 



A \ s 
1 / 



' r(i-c) 

2_ / ' 
A 



.EL 

u 



1 3 
^ + Ye 
A 2 -2 



+ 



A 



.E 

u 



1 

A 



A 

(A 2 -l)F 2 (s,t,u) 



(M^\H^(e)\M^) = {M<®\M<®) 



3^7 



X 



2eT(l-e) 



H (2) 



.EL 

u 



2e 



(4.15) 



(4.16) 



where the square bracket in Eq. (4.16) is a guess simply motivated by summing 
over the antennae present in the quark-quark scattering process and on dimensional 
grounds. Different choices only affect the finite remainder. 
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The functions T\ and JF 2 appearing in Eq. (4.15) are finite functions and are 
obtained from projection of 1^ onto the one-loop amplitude. We find 

Fi{s,t, u) = 
1 



2N 

N 2 -l 



2e(3 - 2e) 



3 «7 



AT 



(N 2 -2)f(s,t,u) + 2f(s,u,t) 
(6 - 7e - 2e 2 ) - (lOe 2 - 4e 3 ) 



1 3 
r(l-e) 2^ 



AT 1 



1 3(1 -e) 



e (3 - 2e) 
iV 2 -l 



s / ^ 
<A^ (0) |A^ (0) > 



/f2 



Bub(s)(7W (0) |7W (0) 
iV 2_2 / ..2 



TV 



t 



2N 



Bub(s) 

[f(s,t,u) - f(s,u,t)\ 



(4.17) 



1 3 
r(l-e) \? + 2~e 



N 



.E. 

u 



N 



t 



(M i0) \M^) 



(4.18) 



where the function f(s,t,u) is written in terms of the one-loop box graph in D 
6 — 2e and the one-loop bubble graph in D = 4 — 2e 



f(s,t,u) 



A(u 2 + t 2 ) - 2e(3ut + Qt 2 + 5m 2 ) - e 2 s(7t + 5m) 



, Qt 2 + 2u 2 - 3es 2 „ B/ . 
+u(l-2e) = Box 6 (s,f). 



Bub(s) - Bub(t) 



(4.19) 



These expressions are valid in all kinematic regions. However, to evaluate the pole 
structure in a particular region, they must be expanded as a series in e. We note 
that in Eq. (4.12), these functions are multiplied by poles in e and must therefore be 
expanded through to 0(e 2 ). In the physical region u < 0, t < 0, Box 6 (w, t) has no 
imaginary part and is given by [17] 



Box 6 (M) = 



3* 2 Ll 



e ^r(i + 6)r(i-6) 2 U 2 

2sY (1 - 2e) (1 - 2e) \ s 
+2e ^Li 3 (a;) - L x Li 2 (:r) 

-2e 2 ^LU(x) + L y Li 3 (x) - iL 2 Li 2 (x) - ±L£ - hlL y + \l\L\ 

+ O (e 3 ) , (4.20) 



7T 2 r2 7T 2 7T 4 \ , 

"T x ~~ ~ 45 J + ( u 



where x = — i/s, = log(x) and L y = log(l — x) and the polylogarithms Li n (z) are 
defined by 



^Li n _i(i) for n = 2, 3, 4 
o t 



Li n (z) 

U 2 (z) = - -log(l-t) 

JO t 



(4.21) 
(4.22) 
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Analytic continuation to other kinematic regions is obtained using the inversion 
formulae for the arguments of the polylogarithms (see for example [22]) when x > 1 

U 2 (x + iO) = -Li 2 Q - i \og 2 (x) + y + m log(x) 

(1 \ 1 7T 2 27T 

-J - - log 3 (x) + y log(x) + y log 2 (a;) 

Li 4 (x + iO) = -LU Q - ^ log 4 (x) + y log 2 (a;) + ^ + y log 3 ^). (4.23) 
Finally, the one-loop bubble integral in D = 4 — 2e dimensions is given by 

Our explicit Feynman diagram reproduces the anticipated pole structure exactly 
and provides a very stringent check on the calculation. We therefore construct the 
finite remainder by subtracting Eq. (4.12) from the full result. 

4.2 Finite contributions 

In this subsection, we give explicit expressions for the finite two-loop contribution to 
^4 8 , Finite, which is given by 

Finite = 2Re(M (0) \M (2)Rn ). (4.25) 

For high energy hadron-hadron collisions, we probe all parton-parton scattering pro- 
cesses simultaneously. We therefore need to be able to evaluate the finite parts in 
the s-, t- and -u-channels corresponding to the processes 

q + q -> 4 + q' 
q + q' -> q' + q 
q + q' -> q + q', 

respectively. In principle, the analytic expressions for different channels are related 
by crossing symmetry. However, the Xbox has cuts in all three channels yielding 
complex parts in all physical regions. The analytic continuation is therefore rather 
involved and prone to error. We therefore choose to give expressions describing 
A 8 (s,t,u), A 8 (t,s,u) and A 8 (u, t, s) which are directly valid in the physical region, 
s > and u, t < 0, and are given in terms of logarithms and polylogarithms that 
have no imaginary parts. 

In general the expansions of the two-loop master integrals [12, 13, 14, 15, 22, 24, 
25] contain the generalised polylogarithms of Nielsen 

W = j^/ ^(')W-0 , „,„>!,,<! (4.26) 

(n — ly.pl Jo t 
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where the level is n + p. Keeping terms up to O (e) corresponds to probing level 

4 so that only polylogarithms with n + p < 4 occur. For p = 1 we find the usual 
polylogarithms 

S n _ 1:1 (z)=U n (z). (4.27) 

A basis set of 6 polylogarithms (one with n + p = 2, two with n + p = 3 and 
three with n + p = 4 is sufficient to describe a function of level 4. At level 4, we 
choose to eliminate the S 2 2, S± 3 and S i2 functions using the standard polylogarithm 
identities [30] and retain the polylogarithms with arguments x, 1 — x and (x — l)/x, 
where 

t u u x — 1 , A 

x = —, y = — = l-x, -- = . 4.28 

s s t x 

For convenience, we also introduce the following logarithms 

L x = log , L y = log , L s = log ^ (4.29) 

where fi is the renormalisation scale. The common choice /i 2 = s corresponds to 
setting L s = 0. 

For each channel, we choose to present our results by grouping terms according 
to the power of the number of colours N and the number of light quarks Np so that 
in channel c 

Finite c = 2(V - lj [n 2 A c + B c + j^C c + N N F D C + ^E c + N 2 F c )j . (4.30) 
4.2.1 The s-channel process qq — > q'q' 

We first give expressions for the s-channel annihilation process, qq — > g'g'. We find 
that 
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We can check some of these results by comparing with the analytic expressions 
presented in Ref. [17] for the QED process e + e~ — > fi + fi~. Taking the QED limit 
corresponds to setting Ca = 0, Cp = 1, T R = 1 as well as setting the cubic Casimir 
C 3 = (iV 2 — l)(iV 2 — 2)/iV 2 = 0. This means that we can directly compare E s ((x 
CfTrNf) and F s (oc T^Np) but not C s which receives contributions from both C 3 
and Cp. We see that (4.35) and (4.36) agree with Eqs. (2.38) and (2.39) of [17] 
respectively. 

The other coefficients, A s , B s , C s and D s are new results. 



4.2.2 The t-channel process q + q' q + q' 

The t-channel process, q + q' — > q + q' is fixed by .4 s (t, s, w). We find that the finite 
two-loop contribution in the t-channel is given by Eq. (4.30) with 
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4.2.3 The w-channel process q + q' — > g + g' 

The w-channel process, g + g' — > g + g' is determined by A 8 (u, t, s). We find that the 
finite two-loop contribution in the w-channel is given by Eq. (4.30) with 
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5. Summary 



In this paper we presented the two-loop QCD corrections to the scattering of two 
distinct massless quarks. Throughout, we have used conventional dimensional regu- 
larisation and the MS scheme to compute the interference of the tree and two-loop 
graphs summed over spins and colours. The pole structure is given in Eq. (4.12) 
while expressions for the finite parts are given for each of the s-, t- and w-channels 
in Sees. 4.2.1, 4.2.2 and 4.2.3 respectively. 

The leading infrared singularity is 0(l/e 4 ) and it is a very strong check on 
the reliability of our calculation that the pole structure obtained by computing the 
Feynman diagrams agrees with that anticipated by Catani through to 0(l/e). For 
the finite N F /N and Np contributions, we agree with prior QED calculations. 

These results form a crucial part of the next-to-next-to-leading order predictions 
for jet cross sections in hadron-hadron collisions. However, they are only a part of 
the whole and must be combined with the tree-level 2 — > 4, the one-loop 2 — > 3 as 
well as the square of the one-loop 2-^2 processes to yield physical cross sections. 
For the most part, the matrix elements themselves are available in the literature. 
Each of the contributions is divergent in the infrared limit and a systematic proce- 
dure for analytically canceling the infrared divergences needs to be established for 
semi-inclusive jet cross sections. Here again, some progress has been made by ex- 
amining the limits of tree-level matrix elements when two particles are unresolved 
[31, 32] and the soft and collinear limits of one-loop amplitudes [33, 34]. The some- 
what simpler case of e + e _ — > photon + jet at next-to-leading order which involves 
the triple collinear limit of tree-level matrix elements as well as the collinear limit of 
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one-loop amplitudes has already been studied [35] and indicates that the technical 
problems are not insurmountable. Another technical difficulty will be to isolate the 
initial state singularities and correctly absorb them into the parton density functions 
at next-to-next-to-leading order. Recent progress towards the three-loop splitting 
functions [36, 37, 38] together with accurate parameterisations in x-space [39, 40] 
suggest that the factorisation can be achieved. Modifications to the global fits to 
provide parton density functions appropriate for next-to-next-to-leading order cal- 
culations are already underway [41]. In summary, we are therefore confident that 
these problems will soon be overcome thereby enabling the analytic cancellation 
of the infrared divergences and the construction of numerical programs to provide 
next-to-next-to-leading order QCD estimates of observable scattering cross sections. 
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